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Abstract 

It has been stated that for a short-ranged surface interaction, the probability of a low-energy 
particle sticking to a surface always vanishes as s ~ k with k — > where k = y/E. Deviations 
from this so-called universal threshold law are derived using a linear model of particle-surface 
scattering. The Fredholm theory of integral equations is used to find the global conditions necessary 
for a convergent solution. The exceptional case of a zero-energy resonance is considered in detail. 
Anomalous threshold laws, where s ~ k 1+a ,a > as k — ► 0, are shown to arise from a soft gap 
in the weighted density of states of excitations; a is determined by the behavior of the weighted 
density of states near the binding energy. 

PACS numbers: 68.43.Mn, 03.65.Nk, 68.49.Bc, 34.50.Dy 
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It has been previously established, both theoretically [l|, |2j, |3j, ]% la, la [7| and experimentally 
that at sufficiently low incident energy E, the probability of a neutral particle sticking 
;o an insulating substrate vanishes as s ~ k with k — > (k = \fE). It has been asserted 
^| that this threshold law is independent of the details of the theoretical model and is thus 
universal. 

It is demonstrated that using a generalization of a linear model of particle-surface scat- 
tering previously considered anomalous threshold laws where s ~ k 1+a , a > as k — > 
are possible. These anomalous laws arise with the presence of a soft gap in the weighted 
density of states of excitations at the binding energy of the particle in the static surface 
potential. Thus, the many-body details of the model can indeed alter the threshold law. 

The model is analyzed using the Fredholm theory of integral equations. There are two 
advantages of this approach: (1) rigorous conditions on the validity of the solution can be 
obtained, and (2) the solution is known to be convergent for all coupling strengths. A gen- 
eral expression for the sticking probability is obtained in terms of the Fredholm determinant. 
Fredholm theory has been previously used in quantum scattering from a static central po- 
tential by Jost and Pais [10j . Here, its utility is illustrated for the case of inelastic surface 
scattering. 

The Hamiltonian is taken to be 



H = Tip + l~C e + T~Cj (1) 

where 

n r = £ + vw, (2) 

H e = ^Uidldi, (3) 

i 

H i = -lJ2 Wi F ( z ) ( 4 ) 

i 

where F(z) = V'(z), and Ti p is the Hamiltonian for the particle moving in the static po- 
tential, V(z). Ti e is the Hamiltonian for the surface excitations; Hi contains the particle- 
excitation coupling; m is the particle mass; Ui is the frequency of the excitation in the ith 
mode; W{ is a generalized mode-dependent coupling; and a\ and a* are excitation creation 
and annihilation operators respectively. 

The Hilbert space is truncated to include only single excitation states. Thus the statistics 
of the excitation does not matter here. The excitation states can be labeled by \n) where n 



represents the mode that is occupied. The system wavefunction can be expanded as 

The following coupled set of equations result 

(H p -k 2 )\4> ) = 7 £m w |&) (6) 

i 

(Hp-e + LuM) = 7M jO |0o) (7) 

where 

M 0i = Wi F (8) 

where F = V and units are chosen so that E = k 2 . 

For the Green's functions to exist for real positive energies, it is neccessary to continue 
the excitation energies into the lower-half of the complex plane uij — » uj — iS. This negative 
imaginary part of the excitation energies corresponds physically to a reciprocal decay time 
that is always present in experiment. It may be thought of as arising from anharmonicity, 
disorder, or interactions with additional degrees of freedom not explicitly contained in the 
model. 

The coupled set of equations can now be converted into a set of coupled integral equations 
with the introduction of the outgoing Green's function Gi = {k 2 — 7i p — oji + if]) -1 . 

/•OO 

(f>i(z;k) = ijjo(z; k)8 i0 + 7^ / dz' K ij (z,z';k)<f) j (z';k) (9) 

3 J ° 

with the kernel Kij(z, z'\ k) = —Gi(z, z'\ k)Mij(z'). The incident wave ipo is placed in the 
zero-excitation (or elastic) channel, a boundary condition consistent with a zero-temperature 
substrate. Normalization of ipo is chosen to give unit flux for the incoming piece. 

Such a system of integral equations can be formally solved using Fredholm theory, pro- 
vided the kernel is sufficiently well-behaved. The coordinate dependence of M is local, and 
the kernel can be modified using a transformation for polar kernels jll| : 

K^z'-k) = -\F(z)\^G i (z 1 z , ;k)M ij (z , )\F{z , )\-^ (10) 

It is straightforward to show that if the following conditions are satisfied, such a kernel 
is in the Hilbert-Schmidt class, 



dz z n \V{z)\ < oo, n = 0,l,2 (11a) 



/>oo 

/ dz z n \F(z)\ < oo, n = 0,1,2 
Jo 



(lib) 



oo i 

^^h'-u + is 



< oo 



(11c) 



where the weighted density of states of excitations is given by p{oS) = \W {uj)\ 2 p{uj) , p{ui) 
being the density of excitation states. The continuum limit has been taken in Eq. Illcl 
These conditions insure that Tr KK* < oo. 

The solution of the channel wavefunctions is expressed in terms of the Fredholm resolvent 
lZij(z,z';k) = r yV i j(z,z';k,j)/A(k,j) of the kernel K. 



(j)i(z) = ip (z)5 i0 + 
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A(* )7 ) 



dz'\F(z)\~ l iV l0 (z, z'; k,j)\F(z')\ho( z '; k) (12) 



Both the numerator of the resolvent and the Fredholm determinant A can be ex- 
pressed as power series expansions in the coupling strength 7. It is well-known that these 
series expansions are absolutely convergent. Thus the Fredholm method is well suited for 
investigating the strong coupling regime. 



Vij(z,z') = Kij{z,z!) 



-7)" 



n=l j r , 



nl 



dz\ 



A = 1 



+ EE 



-7)" 



n=l j„ 



nl 



dz\ 



dz r 



Kij{z,z') K ih (z,zi) 
Kfoj (zi , z) Kj xix (zi,zi) 

Kj n j( z n, z ) 



Kjij n {Z, z n) 



■Kj n j n \ z ni z n) 



dz r 



(13) 



■Kjnji(. Z m z l) Kj n j 2 (z n , Z<i) " " " K ,j n j n ( y Z n , Z n 

(14) 



Since M has the form of a border matrix, these series expansions are truncated. This 
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greatly simplifies Eqs. and El to the following 

V i0 (z,z';k,>y) = K i0 (z, z'\ k) + 7<5 i0 ^ / dziK^z, Z\\ k)K m Q(zi, z'\ k) (15) 

m 

A(fc,7) = l-7 2 ^ / dzdz' K 0m (z,z';k)K m0 (z',z;k) (16) 

m 

Elements of the S-matrix can be obtained from the asymptotic form of the exact solution. 

_ /(*) A(-fc, 7 ) h7 , 

So ° ~ /(-*) A(*, 7 ) (17a) 

2i7v / M^ /" u(/c n , 2) 2;) , , ( . 

s »" = " a(W 7 * 7Pw 7R) (17b) 



where /(&) is the Jost function associated with the potential V, k n = y/k 2 — u n + and 
u(/c, z) is the regular solution in the potential V with the condition that u'(k, 0) = 1. 

In the absence of sticking, the open-channel S-matrix is unitary. Conservation of particle 
flux implies that the deviation from unitarity is the sticking probability s(k, 7). Thus 

3^,^ = 1- J2\S0n\ 2 (18) 

kl>0 

The Fredholm determinant is an analytic function of k in the upper-half of the complex 
/c-plane. It is important to note k = is an exceptional point because of the possibility of 
a zero-energy resonance in the geometry of surface scattering. (A centrifugal potential for 
£ 7^ 0, present in scattering from a point source, is absent in the case of surface scattering.) A 
zero-energy resonance in the potential V permits the Jost function to vanish at zero energy, 
/(0) = 0, leading to a simple pole in the Green's function in the elastic channel. Poles in 
the inelastic Green's functions G m have been pushed off the real axis into the lower-half of 
the complex fc-plane, with the addition of an imaginary part to the excitation frequencies. 

The Fredholm determinant can be written as 

A(fc, 7 ) = l + 7 2 r(A;) (19) 



where 



and 



r{k) = - J dzdz'F(z)F(z')G (z, z'; k)l{z, z'\ k) (20) 

poo 

X(z,z';k)= / dujp(uj)G (z, z'\ Vk 2 -uj + i5) (21) 
Jo 



Substitution of Eqs. I17al I17bl and [THJ into Eq. [TB] yields 

7 2 
|A(*, 7 )| 



s(kn) = IA ,7* , 2 Ref2(r c (fc)-r c (-fc))+ 7 2 (|r(fc)| 2 - |r(-fc)| 2 ) ) (22) 



where 



(k) = - dzdz'F(z)F(z')G (z,z';k) dup(u)G (z, z'; Vk 2 -uo + id) (23) 



'fc 2 

The low-energy regime is now considered. The Green's function Gq(z, z'; k) is analytic in 
the first quadrant of the complex k-plane and can be expanded as a Laurent series about 
k = 0. The series terminates at a simple pole to allow for the possibility of a zero-energy 



resonance 



IJ. Thus, 



G (z, z'- k) = 9 -^l + go ( Zj z >) + gx ^ z ') k + . . . (24 ) 

The discussion now focuses on the case where there is a single bound state in the static 
potential V at energy — k 2 with eigenstate real and normalized. With p suitably smooth, 

lml(z, z'; k) = l^z, z'; k) = -np(k 2 + K 2 b )M z )M z '), ( 25 ) 
Since g± is pure imaginary, for the case where /(0) ^ (g_i = 0), 

Re(r c (A;) - r c (-Jfe)) = 2-nk~p{k 2 + K 2 b )A u k^O (26) 

where 

Ai = J dzdz'F{z)F{z')§ h (z)§ b {z')\m gi {z,z') (27) 

Similarly, 

r(k)r*(k) - r(-k)T*(-k) = ^kpik 2 + k 2 )£i, k -> (28) 

where 

B i = J dz 1 dz[dz 2 dz' 2 F(z 1 )F(z[)F(z 2 )F( y z' 2 )g ( y z 1 , z[)Im gi(z 2 , z' 2 ) 

x ($ 6 ( 2 ; 2 )<l> 6 (4)J r (z 1 , 4) - ^(^^(^T^, 4)) (29) 

and Z r (;z, 2;') = P J °° dup(uj)Go(z, z'\ iy/u). 

Thus, the low-energy sticking probability is of the form 
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where 

C = -2 J dzdz'F(z)F(z')g {z,z')I r {z,z') (31) 

D = J dz 1 dz' 1 dz 2 dz' 2 F(z 1 )F(z' 1 )F(z 2 )F(z' 2 )g (z 1 , z[)g {z 2 , z' 2 ) 

x (l r ( Zl , z[)l r {z 2 , 4) + vr 2 p 2 (^)$ fc (^)$ 6 (^)$ 6 ( 22 )$ 6 (4)) (32) 

Under typical circumstances, p(k 2 + k 2 ) approaches a finite and non-zero limit as k — > 0. 
Then the "universal" sticking law is recovered from Eq. EH for the case of a short-ranged 
potential, with s ~ k as k — > 0. 

Anomalous threshold laws, where s ~ k 1+a (a > 0), are theoretically possible for cases 
where the weighted density of states has a soft gap approaching zero at the binding energy. 
It is noted that a soft gap in the density of states frequently occurs in the vicinity of a phase 
transition. Scattering from a substrate undergoing a at low temperature structural phase 
transition might reveal such anomalous laws. The temperature of the substrate could be 
used to tune the softening phonon mode. In the vicinity of a minimum in the dispersion, a 
soft gap in the density of states appears, with p{oj) oc y/u — uiq in the neighborhood above 
the softened frequency ujq. For the case of a smooth weighting function W, non-vanishing 
at nl, a = 1 and s ~ k 2 as k — > 0. 

The phonon and roton contributions to sticking have been independently measured 
Q] for helium scattering from superfluid 4 He. The roton density of states of superfluid 4 He 
vanishes as the energy approaches the minimum energy to create a roton. It is known 
that the static potential can be substantially modified by preparing a film of superfluid 
4 He on a substrate with a van der Waals coefficient different from that of the 4 He film. 
If the depth of the binding energy is adjusted to the roton minimum energy by using the 
appropriate substrate and film thickness, then it might be possible to measure an anomalous 
threshold law for sticking by roton creation. 

Another possibility would be to use a substrate where the optical phonon energy at the 
zone edge is equal to the binding energy. While surface preparation could be used to adjust 
the binding energy, the van Hove singularity of a bulk optical phonon at the zone edge 
is insensitive to surface perturbations and can provide the necessary soft gap. Flatte and 
Kohn have discussed the general conditions for dominance of inelastic scattering by bulk 
substrate parameters. 
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Consider the case of p(u) ~ (u — kI) v for uj > k\ (77 > 0). The sticking probability then 
behaves as 

The suppression of the weighted density of states at the binding energy results in a decline 
in the sticking probability at low energies that is more rapid than would be predicted from 
the "universal" threshold law. Under the conditions considered, the frequency dependence 
of the weighted density of states near the binding energy appears in the energy dependence 
of the sticking probability near zero energy. 

While a zero-energy resonance is not a true bound state, it may be present in addition 
to a bound state. The presence of the zero-energy resonance affects the energy dependence 
of the elastic Green's function and alters the form of the sticking probability with coupling 
strength 7. 

With /(0) = (<7_i 7^ 0) and 7 7^ 0, at suitably low energies, 

S (k, 1 ) = ^ A -\ + / B ^ k^, k^O (34) 
7 u_i 

where A-\ and are obtained from Eqs. E3 an d I2H by replacing g\ — > g_i and D-\ by 
replacing g —>■ Im g_i in Eq. E21 The same k-dependence of the threshold laws results as 
in Eq. ESI however, the dependence on coupling strength differs. 

Anomalous threshold laws in low energy particle-surface scattering have been shown to 
result from a soft gap in the weighted density of states at the frequency equal to binding 
energy. This provides a counterexample to the claim [9] that the "universal" threshold law 
is independent of model details. Fredholm theory has been used to find an exact solution 
together with global conditions on the potential and inelastic coupling functions for the 
validity of the solution. The solution of the model is valid for arbitrarily large coupling 
strengths 7. A general expression for the sticking probability in terms of the Fredholm 
determinant has been derived. The exceptional case of a zero-energy resonance is seen to 
affect the coupling constant dependence, but not affect the energy dependence of the sticking 
probability. 
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